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^ . Abstract 

The notion of Loschmidt echo (also called "quantum fidelity") has been 
introduced in order to study the (in)-stability of the quantum dynamics under 
^ ' perturbations of the Hamiltonian. It has been extensively studied in the 

past few years in the physics literature, in connection with the problems of 
"quantum chaos" , quantum computation and decoherence. 
In this paper, we study this quantity semiclassically (as h ^ 0), taking as 
reference quantum states the usual coherent states. The latter are known 
to be well adapted to a semiclassical analysis, in particular with respect to 
semiclassical estimates of their time evolution. For times not larger than the 
so-called "Ehrenfest time" C\ log h\ , we are able to estimate semiclassically 
the Loschmidt Echo as a function of t (time), h (Planck constant), and 6 (the 
size of the perturbation). The way two classical trajectories merging from the 
same point in classical phase-space, fly apart or come close together along the 
evolutions governed by the perturbed and unperturbed Hamiltonians play a 
major role in this estimate. 

We also give estimates of the "return probability" (again on reference states 
being the coherent states) by the same method, as a function of t and h. 
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1 Introduction 



The semiclassical time behaviour of quantum wavepackets has been the subject 
of intense interest in the last decades, in particular in situations where there is some 
hyperbolicity in the corresponding classical dynamics (Lyapunov exponents) [9], [T7] . 
|30j . Moreover the response of a quantum system to an external perturbation when 
the size 6 of the perturbation increases can manifest intriguing properties such as 
recurrences or decay in time of the so-called Loschmidt Echo (or "quantum fidelity" ) 
[7j, [8]. By Loschmidt Echo we mean the following: 

starting from a quantum Hamiltonian H in L^(R"'), whose classical counterpart H 
has a chaotic dynamics, and adding to it a " perturbation" Hs = H + SV, then 
we compare the evolutions in time U{t) := e~**^/^ , Us{t) := e"**-'^*/^ of initial 
quantum wavepackets (p sufficiently well localized around some point z in phase- 
space; more precisely the overlap between the two evolutions, or rather its square 
absolute value, is: 

For example for quantum dynamics in Hilbert space 7i = L^(R'^), d being the 
space dimension, ip can be chosen as the usual coherent states, since they are the 
quantum wavepackets "as most localized as possible" in phase-space M^'^. 

Since for 5 = 0, we obviously have FfiQ{t) = 1, and for any 6 ,Ffis{0) = 1, 
the type of decay in t of Ff^^t) so to say measures the (in) fidelity of the quantum 
evolution with respect to a perturbation of size 5 for generic initial wavepackets p. 

The notion of Loschmidt Echo seems to have been first introduced by Peres 
( |24j ) ■ in the following spirit: since the sensitivity to initial data which characterizes 
classical chaos has no quantum counterpart because of unitarity of the quantum 
evolution, at least the "sensitivity to perturbations" of the Hamiltonian could replace 
it as a characterization of chaoticity in the "quantum world" . 

A big amount of recent work appeared on the subject, studying in an essen- 
tially heuristic way the decay in time of Ff^^t) as t increases from zero to infinity; 
some of them also study this point in relationship with the important question of 
decoherence. (See P, [5], [12-14], [18], [22-26], [30-32]). 
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In this "jungle" of sometimes contradictory results, it is hard to see the various 
arguments involved, in particular the precise behaviour of Ff^sif) as 5 (the size of 
the perturbation), t (the time), and of course h (the Planck constant) are varied, 
in particular in which sense and order the various limts 5-^0, /I— >0,t-^cx3 are 
taken. 

Also an important point to consider is how Fji^^{t) depends on the location of the 
phase-space point z around which the initial wavepacket if is peaked (since classical 
chaoticity distinguishes various zones in phase-space with "more or less regularity 
properties"). 

The aim of the present paper is to start a rigorous approach of the question of 
semiclassical estimate of Fn^s{t)i in terms of classical characteristics of the (perturbed 
and unperturbed classical flows), for initial wavepackets = being the coherent 
state at phase-space point z. These estimates are non-perturbative, and are carefully 
calculated in terms of parameters {z,6,t,h). The main tools we have used and 
developed in this respect are 

1) semiclassical coherent states propagation estimates ([9]) 

2) a beautiful formula inspired by B. Mehlig and M. Wilkinson {[221) about the Weyl 
symbol of a metaplectic operator, and thus of its expectation value in coherent states 
as a simple Gaussian phase-space integral ( see [10] where we have completed the 
proof of Mehlig- Wilkinson, and treated in particular the case where the monodromy 
operator has eigenvalue 1). 

Note that very recently, J. Bolte and T. Schwaibold have independently obtained a 
similar result about semiclassical estimates of the Quantum Fidelity ([2])- 

The plan of this paper is as follows. In section 2 we give some preliminaries about 
the Echo for suitable quantum observables, and give the semiclassics of it. In Section 
3, we consider the (integrable) d = 1 case, and consider the "return probability" in 
the semiclassical limit. We give a mathematical rigorous presentation of beautiful 
results on "quantum revivals" obtained by physicists twenty years ago (see [23], [32] . 
|20j). In Section 4 we consider the general o?- dimensional case and give a semiclassical 
calculus of the "return probability" and of the quantum fidelity, with precise error 
estimates. 
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2 Preliminaries 

Let us consider the quantum Hamiltonian Hg = Hq + 6V, depending on a real pa- 
rameter 6. Hs is the Weyl quantization of smooth classical observables defined on 
the phase space M^"'. Our assumptions on are as follows: 
Assumptions : 

(HI) Hs e C^iM?'^) and \dlHs{X)\ < C^, VX e B?^, V7 : I7I > 2 
or 

(H'l) There exist a bounded open set f2 C R^'* such that VL is left invariant by the 
classical flow 0^, Vt, 5 G M defined by the classical Hamiltonian H^. We assume that 
the ^-Weyl quantization = Op^Hs of Hs is a self-adjoint operator in L^(M'^) for 
all 5 eR. 

Let L G C°^(M2"') be a classical observable and L = Op'llL. Then we assume: 
(H2) L G 5(M2'^) if (HI) is satisfied, 

(H'2) L G C^{Vl) if (H'l) is satisfied, where C^{Vt) is the linear space of C°°-smooth 
functions with compact support in fl. 

Let us consider the time evolution unitary operator Us{t), in the Hilbert space 
H = L2(M^), 

/ it . 
Us{t) = exp --Hs 



h 

Definition 2.1 (i) The quantum echo is the unitary operator defined by 

E^'\t) = Uo{-t)Us{t) (2.1) 
(a) The quantum fidelity, for a state ipQ, \\ipo\\ = I, is defined by 



/f(t)= {^o,El'>m,) (2.2) 

(Hi) Let Ls{t) be the following time- dependent quantum observable: 

Ls{t) := Elitr'LElit) 

The notion of "fidelity" was introduced first in classical mechanics by Loschmidt (in 
discussions with Boltzmann) then adapted in quantum mechanics by Peres [24] . 



n(9)/ 
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Let us define 0^, the classical flow defined in the phase space Z = R^'^ by the 
classical Hamiltonian Hs- Recall that zs^t '■= (Psi^o) is the solution of the differential 
equation zt = JVHs^zt), zt=o = zq. So that the "classical echo" is defined by 

Here J is the symplectic matrix given as: 

J :-( ' ) (2.3) 



, -Id ^ 
and Id is the identity d x d matrix. 

We can see easily that in the semiclassical limit, h ^ 0, the quantum echo converges 
to the classical echo. In more mathematical terms, the quantum echo is a ^ - 
Fourier Integral Operator whose canonical relation is the classical echo. This is a 
consequence of the semiclassical Egorov theorem as we shall see now, at least when 
the reference quantum state is a "coherent state". Let us recall here the definition 
of a (Gaussian) coherent state which will be used later: 
Given ipo{x) := (vr/i)""'/^ exp(— we define, for z := {q,p) € M^*^: 



(Pz '■= T{z){pQ, where f{z) := exp 



i{p.Q - q.P) 

h 



are the translation Weyl operators. 

Proposition 2.2 Assume either (H1-H2) or (H'l-H'2) for the Hamiltonians 
and observables L. 
(i) We have 

\im{2TThYi:i{LL5{t)) = [ L{X)L{Ef\t,X))dX 
(a) Let be ipz the coherent state living at z. Then we have: 

hm{El'\t)^z\LEl'\t)vz) = L{Ef\t, z)) (2.4) 

Moreover there exists C > such that the limits are uniform as long as \t\ < 
C\ \ogh\. If Hq is integrable in and if S = 0{h), then the limit is uniform as long 



as 



\t\ < ah-^/^+' 



Proof: It follows from the semiclassical Egorov Theorem,with improvement for 
large times derived by Bouzouina- Robert ([3]). 
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Remark 2.3 An important question is to control the time of validity of the semi- 
classical approximation. Rigorous mathematical results are far from numerical and 
theoretical expected physical results. Without assumptions on classical flows this time 
is the Ehrenfest time (of order log(^~^)). 

3 Revivals for 1-D systems 

In this Section we shall consider the Return Probability, which is a simplified form 
of Quantum Fidelity as we shall explain in the next Section. For one-dimensional 
problems, this Return Probability manifests interesting recurrences very close to 1, 
as time evolves. This phenomenon was studied in the physics literature to under- 
stand time evolution of Rydberg atoms and their quantum beats, with decay and 
reformation of the wave packet (see for example [231 120]. [32] and references herein 
contained). In this section we want to give a fiavour of results obtained by physicists 
in the last twenty years, concerning revivals for the quantum return probability ([32] 
for a very clear and detailed review) and show how to put them in a more rigorous 
mathematical framework. Let us consider a classical 1-D Hamiltonian H. One as- 
sumes if to be a smooth, confining with one well Hamiltonian. This means that the 
energy surface H~^{E) has only one connected component in phase-space Z. Let 

be an orthonormal basis of eigenstates, with eigenvalues En, n & N. 
Let ipo = ^^c„\E'„, an initial normalized state, and ipt = U{t)ipQ. Then the autocor- 

relation fonction is : 

a(t):=(^o|^t)=$^|cn|Vf^" (3.1) 

neN 

and the return probability is defined by 

p{t) = \a{t)\\ (3.2) 

Let us remark here that a is an almost periodic function (in the sense of H. Bohr) in 
time t on M. Therefore, for every e > 0, there exists > and for every A; G Z there 
exists tk G [A;T£, {k + '^)T^[ such that \a{tk) — 1| < e. This fact can be interpreted as 
a quantum analog of the famous return Poincare in classical mechanics. 
But we have no information here on the almost return time t^, in particular when h 
tends to zero. For 1-D systems much more accurate results are available because for 
these systems the spectrum can be computed with error 0{ff^) according the Bohr- 
Sommerfeld quantization rule. Recall here this result. We take the presentation 
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from the paper by Helffer- Robert ([H]) and we refer to this paper for more details, 
(see also the thesis of Bily for a proof using coherent states). 
Let us give now the sufficient assumptions: 
[Ai) H{z) IS real valued H e C°°(M2). 

{A2) H is bounded belovqj : there exist Cq > and 70 G M such that Cq < H{z) +70. 
Furthermore H{z) + 70 is supposed to be a temperate weight, i.e there exist C > 
0, M eR, such that : 

H{z) + Ao < C{H{z') + Ao)(l + \z- z'\)^ Wz, z' G Z. 

{A3) V7 multiindex 3c > such that: \d2H\ < c{H + Aq). 

We want to consider here bound states of -ff in a fixed energy band. So, let us 
consider a classical energy interval / =]E_ — e,E^ + £[, E_ < such that we 
have: 

{A^ H^^{I) is a bounded set of the phase space M^. 

This implies that in the closed interval / = [E_,E^], for h> small enough, the 

spectrum of in / is purely discrete ([T8]). 

For some energy level E E I, let us introduce the assumption : 

(As) E is a regular value of H. That means: H{x,C,) = E ^ V[x,^)H{x,^) 7^ 0. 

Furthermore we assume that for every E E I, H~^{E) is a connected curve. 

Let us consider a non critical energy interval [E^,E^]. It is well konwn that the 

action integral is Jl{E) = /f^^^-j^^ dz and the period along the energy curve H~^{E) 

is Te = J'{E), E G [E~,E+]. Let us denote F± = J{E^). The eigenvalues of H in 

[E^,E'^] are determined by the following Bohr-Sommerfeld rule. 

Theorem 3.1 (| |18] ) Under the assumptions (Ai) to (A^), there exists a C°° func- 
tion on [F~,F~^], F \—>- b{F,h) and functions bj defined on [F~,F~^] such that 

b{F, h) = ''^^bj{F)h^ + 0{h^) and the eigenvalues En of H in I are given by 

jen 

En = b{{n + -)h, h) + C»(n~), for n such that {n + -)h G [F", F+] (3.3) 
2 2 



where 60(F) = 2'kJ-^{F), and 61 = . (3.4) 



Remark 3.2 In recent papers JD^, JH]/ the authors have given some methods to com- 
pute explicitly the terms bj for j >2 in the expansion in h in the Bohr-Sommerfeld 
rule. 



^ Using the semi-classical functional calculus |29] it is not a serious restriction 
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Let us now choose an initial wave packet, = ^^c„\l/„, tightly spread around the 

n 

energy En where En = b{n + ^,h) and n is a family of given quantum numbers, 
depending on h, and such that (n + G for every h g]0, /iq]- Let us 

choose the coefficient Cn defined by 

f En — En\ f Er, — E 

Cn = Kr^hXl XO 



where xi ^ '^(M), xo has a bounded support, supp[xo] ^] — ^,^[, Xo(a^) = 1 
on [—1/2,1/2], and K^-fi is defined such that the L^-norm of the wave packet is 
|c„p = 1. We shall choose Tf^ and e^^ such that ^ = 0{h^), for some 5 > 0. 

neN 

Remark 3.3 From a physical point of view, a state ip above is prepared by ex- 
citing an atom with a laser beam. The new object is a Rydberg atom. 

Pratically, we shall choose r^i = h^. We define: cr = ^ = hP^-^ and = h^' , with 
< 9' < 9 < 1. It is more suggestive for us to keep the notations r^^ abar and e^,. 
Let us apply the Taylor formula to b{F, h) around F = {n+ ^)h. 

En -En = hb'oin -n) + ^b'^{n - n)^ + ^b'^{n - n)^ + h%{n - n) + O(r^), (3.5) 

Z 

if In — TT-I < Ca, with C > and where the derivatives of bj in F are computed in 
F. Up to a small error in h, it is possible to change the definitions of Xi s'^nd K^-^fi 
such that Cn = Kr^nX ("^v^)) with a = ^. Let us remark that Kt-^ is of order a^^. 
This is easily seen from the following lemma. 



Lemma 3.4 With the previous notations and assumptions we have 

\Xi{x)\'^dx (3.6) 




E„ — En\ f En — E^ '"^ 



n 

X\ I I Xo 



Proof. Besides the assumptions, we make use of formula (13.51) and of the following 
well known estimate for the number of bound states 



#K Ene[E-e,E + e]} = 0(^^) 



The details are left to the reader. □ 

Let us denote by ai{t) the approximation for a{t) obtained by plugging in (13. ip 
the i-first terms of the Taylor expansion (13.51) denoted by Ki{n) (1 < i < 3). So we 
get the following preliminary result: 
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Proposition 3.5 We have 

|a(t)|2 = |a,(t)|2 + C(|t|rV'+^) (3.7) 

In particular, \ai{t)\^ is a semiclassical approximation for \a{t)\^ valid for times t 
such that \t\ is less than h}'^^T~^~'^ for any £ > 0, with a reminder term 0{h^). 

Corollciry 3.6 For every e > 0, we can choose < 1, close enough to 1, such that 

\a{t)\'' = \a,{t)\^ + 0{\t\K-') (3.8) 

Prom the proposition and its corollary we can give a mathematical proof for the 

collapses and revivals phenomenon concerning the return probability p{t). 

Let us remark first that Ki{n) — h'Q{n — n) so |ai(i)p is periodic with period = |r 

(classical period along the orbit of energy E). So the return probabihty p{t) is close 

to 1 for i = NTa as far as \t\ is less than r-^h^+^ = 0{h-^+'') {e' > e.) 

For larger times, we have to consider K2{n) — bQ{n — n) + |6o(^ — nY and a second 

time scale dependent on h, the revival time, defined as Trev — Let us introduce 

the integer N — ^^tT\ as a new time origin, lit — NT^ + s we have 

"2(^) = X] |c^+n|^exp (^^^^^ (^~^^ ~ OTce)rnj , 

where 9 e [0, 1[. Therefore we have 

a{NTce + s) = a{s) + 0{W), as long as |s| < h'~'^ (3.9) 

So, around the time NT^, which is of order h^^, the signal retains its initial form 
and moves according to the classical laws. 

Now we shall prove that for large time intervals, below the time Tj^^^ the signal a{t) 
is very small. Let us consider 

where 5i > 0, 82 > are any small fixed real numbers satisfying 52 + Si/2 + 9 < 1. 
Proposition 3.7 Under the previous assumptions and notations, we have 



(3.10) 
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Proof. For simplicity, we shall prove the collapse property ( ]3.10p for a Gaussian 

cut-off, Xi{^) = We can assume that = 27i. 

The trick here is to apply the Poisson formula in the time variable to 

a2(t) = Krfi 2J exp { -— + 2mtj^ \ exp{itm) (3.11) 

mgZ ^ rev/ 

So, applying the classical formula for the Fourier transform of a Gaussian we get 



02 



{t)=K^J—T^xp (-2rr'^-^^) (3.12) 



where ^tA= ^ " 

' \ C' -t rev 

We have 



rev 



and each Gaussian term in the sum in (13.121) has width 6t, given by 

/ 1 f2^2 \ 1/2 

5*=(5?(7,.-^)))-^/^= A + 16vr^^ 

From formula (13.121) . we can see that a sufficient condition for t to be a collapse 
time for p{t) is that and 5t tend to with h. Therefore we get easily (l3.10p .D 

Remark 3.8 The length of Jji is of order is of order h^'^^^"^ . So the length of Jn, is 
very large for h very small (remember that 82 is small and 9 close to 1 ). Therefore in 
the large intervall Jfi, a{t) is very small and in particular its classical period Td has 
disappeared. But we have seen that this period appears again at time NTce (close to 
Trev, for small h. These facts justify the name "revival" given to this phenomenon. 
As it is shown in [32], it is also possible to observe fractional revivals, using some 
elementary properties of integers. 

Remark 3.9 The above analysis could be extended to completely integrable systems 
in d degree of freedom, using the corresponding Bohr-Sommerfeld rules 

Remark 3.10 In the next section, for d-multidimensional sytems, we shall start 
with a Gaussian coherent (pz of classical energy E = H{z). Let us consider x 
above and such that x = 1 m a small neighborhood of E. Then, modulo an error 
term 0{h°°), we have easily 

{^.\U{t)^z) = J2x (^^) \{^.\en)\'e-'^''- (3.13) 
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We get something similar to the definition of a(t) but with coefficients Cn not neces- 
sary smooth in the variable n, so application of the Poisson formula seems difficult. 

4 Fidelity on coherent states 

Let us recall that the Return Probability for suitable time-dependent Hamiltonians 
H{t) (for which the quantum unitary evolution U{t, s) can be shown to exist) in 
some quantum state ipo E H = L^(M'^) is defined as 

R{t) := \{tPo,U{t,0)iJo)\' (4.1) 

It measures the quantum probability that the time-evolved quantum state U{t, 0)ipo 
returns close to its initial quantum configuration i/jq. 

The Quantum Fidelity (I2.2p can be related to the Return Probability for a suitable 
time-dependent Hamiltonian: 

m = exp j (Hs - Ho) exp (^-^ j (4.2) 

Namely, according to Schrodinger equation, we have that, defining U(t) := Uo{—t)Us(t), 

ihjUit) = H{t)U{t) 
so that U{t, s) = U{t)U{—s) is the time evolution associated to (14. 2p . 

Thus as a training for studying Quantum Fidelities, let us first consider the semi- 
classical study of the Return Probability in the coherent states. 

Let us assume the following hypotheses: 



1. H{t,X) is a smooth time dependent Hamiltonian, continuous in time t G M, 
C°° in X G M^'^ such that for every multiindex a there exist Cq, > and 
M„ G M such that 

\d'^H{t,X)\<Ca{l + \X\f^-, forX GM^^tGM. (4.3) 

2. The classical flow 0*''' generated by H{t) (with initial data at s) exists for all 
times t, s. We shall denote 0* = 0*'''. 
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3. H{t) := Op^H(t) is self-adjoint on L'^iM.'^), and generates a strongly con- 
tinuous evolution operator U{t,s) satisfying the chain rule U{t,T)U{T, s) = 

Note that sufficient conditions for this to hold are given in [H]. 

Then we define the Stability Matrix F{t) for the fiow 0* as follows: 

It is the 2d x 2d symplectic matrix solution of the following linear problem: 

Ft = JH"{t)Ft 

where H"{t) is the Hessian of H{t) taken at point (p^z of the phase-space trajectory, 
starting with initial phase-space point z G M^'^, J being the symplectic matrix given 
by (Q. We have: 

Theorem 4.1 Let us assume Hypotheses 1,2,3 above. Then we have, for the am- 
plitude of the return probability r{t,z) := \{U{t,0)ip;^, Lp;,)\, the asymptotic formula 
as h 0, 

r{t, z) = |det(\4)r'^'e^ + 0{Vh) (4.4) 

where 

Vt = ^{l + Ft + tJ{l-Ft)) 
Ft being the stability matrix for the flow, and 

At = h^pX^t - z) ■ {zt - z) 

with Tf, = {1 + iJ){l + Ft){2Vt)-\l - tJ) - II 

In particular if z lies on a periodic orbit 7 of the classical flow, with period T^, and 
if F{T^) is unitary, we get: 

r{T.,,z) = l + 0{h'/^) 

namely we have almost "quantum revival" when h ^ 0. 

The proof will be very similar to the one we establish below for the Quantum 
Fidelity. 

Let us now consider the fidelity problem. We want to analyze 

f!it) :=Kf/o(tV. , Usit)^.)\'. 
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For the generators of the (time independent) Hamiltonians Hs, we assume: 
(HO) Hs = Op-i/5 

(H2) Hs is a smooth Hamihonian such that there exists for any multiindcx 7 
constants > 0, and G M such that \dlHs{X)\ < C^{1 + |X|)™^ VX G M^^. 

Wc denote by Fs{t) the stabihty matrix for Hs (and similarly for Fo{t) and Hq), 
and by 05 (t) the classical flow for Hs, so that the phase-space point of the classical 
trajectory starting from z G M^'' is zf = 4>sZ. 

Then we have: 

Theorem 4.2 Assume (HO), (H2). 

(i) Then for any N > 1 we have the asymptotic expansion: 



where 



ej{t,z,X) := aj{t, z, X)e 



• aj(t, z, .) is a polynomial of degree < 3j 
• At^z '■— -^Fq^TfFq^ , F being the transpose of F 
• Ff := (1 + iJ)(l + F){1 + F + i J(l - F))-^ (1 - i J) - 1 
• and F denotes F :— Fq^Fs 

Rl^\t, h) = 0(/i*^^+-'^)/^) is uniform on every interval [—T,T] (0 < T < 00). 
In particular we have: 

/ n \\-V2 

eo{t,z) = f det f -(l + F + iJ(l-F)) 
and denoting by Vf the following 2d x 2d matrix: 



Vf := ^(l + F + iJ(l-F)) (4.6) 



flit) = I det Vf\-^ exp 0KA,,,(4 - z^) ■ {zf - z°)^ + 0{Vn) 
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(a) Moreover, we have, in the sense of quadratic forms the following inequality: 

1 



-1 r-i 



< ^^^n^'^n 

' 2 + 2||FP ° ° 

where \\F\\ is the largest eigenvalue of F, and for any symplectic matrix F: 



(4.7) 



rmnim/f(t) = l 



I detVFl > 1 
zf = z^ and F = Fq^Fs is a unitary matrix. 



Theorem 4.3 Let us assume that Hs{t) obeys a Gevrey class assumption: 
(G,) : 3p > 0, z/ > 0, C > 0, i? > such that: 

V7 e N'^ \dlHs{X)\ < Ci?^+I^l(7!)^e"l^l'^\ VX e C'^ \^X\ < p, W5 E [0,5o], > 

then with the notations of Theorem 4-2, we have: 

(z) \ej{t,z,X)\ < C^+\j + l)^^^e^^l^l, where = 2s - 1. 

(a) Ve > 3a > 0, 3c > such that for \X\ > e and j < j^prz we have \ej{t, z, X)\ < 

c 

e ft-'-''''* . 



[Ill, 



3C > 0, c> such that if Nn = [j^j^] we have Vt : \t\ < T: 



(f/o(t)^. , \J,{t)^,) - [t, z, j 



< C exp 



The proof of these theorems heavily rehes on a resuh for semiclassical propaga- 
tion of coherent states (see f9]) which has been revisited in [30], [31] : 

Theorem 4.4 Assume (HO), (H2). Then there exists a family of polynomials {hj(t,x)}j(z^^ 
in d real variables x = {xi, ■ ■ ■ , Xd), with time dependent coefficients, such that for 
all h g]0, 1], we have 



U{t)^, - exp ( ^ j f{zt)KnRi{Ft) ( h^%{t)g 

0<j<N 



)/2 



L2 



(4.8) 

such that for every N E N, and every T < +00 we have sup C{N, t, z, h) < +oo. 

0<h<l,\t\<T 

g is the simple normalized gaussian function M: 



g{x) := (7r)-^/^exp(- 



X 
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Here Ri{F) is the usual metaplectic representation (for h= I) associated to F (see 
fTD^ )- Moreover Afi is the following unitary transform in Ti: 

AaV^(x) = h-'^'^ij {xh-^'^) and 



1 



lt{z) = - Zs- VH{zs)ds - tH{z) 
^ Jo 

Let us here recall a simple property of the metaplectic representation: 

if 

\C D 

is the 4 d X d block- matrix form of the symplectic matrix F, the action of Ri{F) 
on the state g is given by: 

R^{F)g = 7i~'^^\det{A + iB))-^/^ exp ■ 

with r := (C + iD){A + iB)~^ 

Let us denote by ip^^'' the approximation of U{t)(pz given by fl4.8p . 
Let us recall some more accurate estimate obtained in [9] and [30 



(i) Let be fixed and R> such that \zt\ < i?, Vt G M. Then there exist cat > 0, 
kR> such that 

h^^+^y^C{N, t, z, h) < CNkR (^^\Ft\^y^^ (1 + \t\f+^ (4.9) 

In particular, in the generic case, we have a positive Lyapunov exponent 7 
such that I -Ft I < e*^'*', so that the semiclassical approximation is valid for 
|t|<^|log^|. 



In the integrable case we have \Ft\ < c\t\ and the semiclassical approximation 
is valid for |t| < h^^^^^^, for any e > 0. 

(ii) If H satisfies the following analyticity assumption in the set 

Qp = {X eC^",\QX\< p} (4.10) 
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where '^X = C^Xi, ■ ■ ■ , '^X2d) and | ■ | is the Euchdean norm in R^'^ for the 
Hermitean norm in C^'^. So we assume there exist p > 0, C > 0, u > 0, such 
that H is holomorphic in Qp and for all X G fip, we have 

< Ce'^l^l. (4.11) 
Then the X-dependent constant in (14.91) can be estimated by 

cjv < C^+'(iV + l)^ (4.12) 

From this estimate we get an approximation for U{t)ipz niodulo an exponen- 
tially small error (see also [T7]). 

(iii) There exist t>0, a>0, k>0 such that for X = {|} (the nearest integer to 
|), we have 

< fce-^ V;ie]0,l]. (4.13) 

Now we apply the above estimates and the results already proven [lOj concerning 
the action of metaplectic transformations on Gaussians. Our aim is to study the 
fidelity 



fsAt) = \{Uo{t)v>z,m)fz)\' (4.14) 

We shall add the index 6 to keep track of the dependence on the perturbation 
parameter in the Hamiltonian Hs- z is fixed so we shall omit index z. 
We use the approximants and ipi'^^Q for both terms of the scalar product in 

fl4.14p . This yields that to get the result of Theorem 4.2 mod 0{Vh) we have to 
calculate: 



f (^=1) (^^) R\F,')g , R\F,%'j (4.15) 
But f l4.15p is simply of the form 



where 

1 

— ( 
2 

Recall that a{X, Y) = X ■ JY is the symplectic form in R^*^. 



P,:=--aizt,zf^ 
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We are just left with the calculation of the matrix elements of the metaplectic 
operator between two different coherent states. We have thus established the follow- 
ing important result [TU] (here everything is independent of h and we have denoted 
gj,:=f(^=^){X)gy. 

Lemma 4.5 The matrix elements of R^{F) on coherent states gx are given by the 
following formula: 

{gy+^\R\F)gy_.) = 

2'^ (det Vt)~^^^exp S^{Kf -1)Y-Y + ^a{X, Y - KpY - KpY) + ^JKpJX ■ x| 

(4.16) 

where Kp := (1 + F){1 + F + iJ{l — F)) and Kp is the transpose of Kp. 

Remark 4.6 //det(l + F) ^ 0, then 

Kp= (]l + a(l-F)(ll + F)"^)"' 

In this case ix{F) has a smooth Weyl symbol given by the following formula: (see 
[TD^ where we have named this formula the "Mehlig- Wilkinson formula", according 
to the physics literature f2^) 

R{F,X) = e'^" \det{l + F)\-^^^ exp {-iJ{l - F)(ll + F)-^X ■ X) (4.17) 

where v is the Maslov index that we have computed exactly. 

Moreover we have: 

Lemma 4.7 For any symplectic matrix F , consider the matrix Vp defined by Ili4.6\ ). 
We have that \ det Vf\ > 1, and \ det Vf\ = 1 if and only if F is unitary. 

Proof of Lemma 14.71 Let 

F = 

be the 4-block decomposition of the 2d x 2d symplectic matrix F. We have the 
following diagonalization property of the Hermitian matrix iJ: 

a = t/(- 

where U is the unitary matrix 



[c d) 




Thus we have: 

Vf = Iu(( M U*FU + f ? ? ]]U* 



2VVooy 'vo2 

and therefore 

= d.t 1 ( ^ + ^ -C) A-D -HB + C) ^ ^ 1 ^ ^ ^ ^^(^ _ 

We conclude that (recall that A is the transpose of the matrix A): 

I det Vpp = det ^[i + D - i{B - C)][A + D + i{B - C)] = det[l + L*L] 

with 

L = ^[A- D + t{B + C)] 
where we have used the symplecticity of F, namely that 

AC-CA = DB-BD = 

AD-DA = 1 



□ 

End of Proof of Theorem 4.2: Putting X = 2Y = ^Fo t{zo,t - zs,t) in KTB . we 
get (US}. 

Now the estimate (14 .Zp easily follows from the following: 
Lemma 4.8 Let 

Then for any X G M^'^ we have: 



2(1 + Si.) 



where sp = is the largest value of FF (F being the transpose of the matrix 

F). 
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Proof: let us begin to assume that det(ll + F) ^ 0. Then we have: 

Kp = (1 + iNy^ where N = J{1- F){1 + is real symmetric 
so we can compute 

^{Kp) = (1 + A^^)"^ = KfK*p and '^{Kp) = -N{1 + N^)'^ 

So we get: 

^-fF{X) = 1 ((1 + JN)KfK*p{1 - NJ)X ■ X - 2|X|2) 
By definition of Kp, we have: 

(1 + JN)Kp = 2 ((1 + iJ)F'^ + 1- ijy^ := 2Tp 
We have, using that F is symplectic 

{T*p)-^Tp^ = 2{F-^F-^ + 1) 

Hence we get: 



TpT; -- = (2{F''F-' + 1) 



2 2(1 + FF) 

IV12 1 ~ 1 

TpTtX ■ X - = — (1 + FFY^X ■ X < -IXP 

^ ^ 2 2^ ^ - 2(l + s^)' ' 

and the conclusion follows for det(ll + F) 7^ 0, hence for every symplectic matrix F 

by continuity. 

□ 
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